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Visible to near-infrared reflectance spectra of macroscopic mixtures have been shown to be linear combinations of 
the reflectances of the pure mineral components in the mixture (1). However, for microscopic or intimate mixtures the 
mixing systematics are in general non-linear (2,3). The sysfcematics may be linearized by conversion of reflectance to 
single-scattering albedo (w), where the equations which relate reflectance to w depend on the method of data collection 
(bihemisphericial, directional hemisphericial, or bi-directional reflectance) (4). Several recently proposed mixing models 
may be used to estimate mineral abundances from the reflectance spectra of intimate mixtures. These models are 
summarized below and a revised model is presented. 

Johnson et al (5) present a semi-empirical mixing model based on the work of Hapke (4). They used directional- 
hemisphericial reflectance data measured in the laboratory and combined pure mineral spectra to fit the spectrum of a 
mixture using an iterative least squares approach. Although the fits of endmember spectra to mixture spectra were 
within experimental error, the predictions of mineral abundances were 2% to 30% from the actual values. A more recent 
approach has been to apply principal components analysis (PCA) to collections of laboratory and lunar spectral data 
measured as or, converted to directional-hemisphericial reflectance (6,7). This technique may be used to determine the 
type of mixing (macroscopic vs. intimate), identify the potential endmembers from a suite of spectra and estimate the 
relative proportions of endmembers in a sample. The advantages of PCA are the identification of endmembers with 
limited a priori knowledge and approximate mixing relations between endmembers. The disadvantages are that it 
requires careful interpretation of die principal axis of variation and the quantitative determination of endmember 
abundances becomes complex for mixtures containing more than two components. 

We have used a similiar curve fitting approach as Johnson et al (5) with the following exceptions. We use a non- 
iterative ( i.e. linear) least squares approach and the data, measured as bi-directional reflectance with incidence and 
emergence angles of 30 c and 0° were converted to w by the following simplified version of Hapke’s (4) equation for bi- 
directional reflectance 

r(i,e)= w (2.2 + H( p )H( M Q )) H( p )= 1 + 2 p (1) 

(i+x 0 ) 1 + 2 M7 

where r is the measured reflectance, i is the incidence angle, e is the emergence angle p =cos(e), p 0 = cos(i) and y 
— (l-w)* 1 . Hapke (4) predicts that for intimate mixtures, the mean single-scattering albedo is related to the single- 
scattering albedo of each of the components by 

w(X)=(|>i(X )M i p i d i )/( l Mpjd;) (2) 

where Mj, pj, and d- are the mass fraction, density and mean particle diameter of each of the i-components in the 
mixture. The endmember mineral spectra were combined to fit the measured mixture spectra by a single least squares 
inversion of equation (2). The fractional contribution of each endmember is given by the relative geometrical cross 
section (f) where ^“(Mj/pjd:)/( T M^pj^) for the jth component of an n -component mixture. 

This model was tested* with two mixture series composed of 45-75 um particles: an anorthite-enstatite series (two 
minerals of similiar overall albedo) and an olivine-magnetite series (two minerals with strongly different overall albedos). 
The quality of the least squares fits as determined by the standard deviation are not greater than 5 k 10“^ in all cases. 
The standard deviation is determined by curve shape only and is not sensitive to the effects of sample preparation which 
may result in up to a 2% offset in the measured reflectance curves (5). Shown in Fig. 1 are the observed and least 
squares fit spectra for the anorthite-enstatite mixture series. 

The computed and actual values for f are shown in Figs. 2 and 3. It is evident from Fig. 2 that, for mixtures of 
comparable albedo, equations (1) and (2) can be used to predict the known fractional abundances to within the 
experimental error. This degree of accuracy is more than sufficient for the purposes of estimating surface mineral 
abundances from known or assumed endmembers. The simplified version of Hapke ’s (4) equation is less satisfactory for 
mixtures containing low albedo materials as shown in Fig. 3. The model consistently overpredicts the olivine content and 
underpredicts the magnetite abundance. The systematic deviation of computed from actual values of f for the mixtures 
containing low albedo minerals suggests that either a more rigorous treatment of Hapke’s (4) equations is needed or an 
empiricial adjustment to equation (1) will allow a more accurate solution to be developed. 

These data indicate that the simplified version of Hapke ’s equation for bi-directional reflectance (equation (1)) may 
be used to deconvolve reflectance spectra into mineral abundances if appropriate endmembers are known or derived from 
other techniques. For surfaces that contain a significant component of very low albedo material, a somewhat modified 
version of this technique will need to be developed. Due to the nature of equation (2) above, this model can also be easily 
expanded for use with multicomponent mixtures. Since the abundances are calculated using a non-iterative approach, 
the application of this method is especially efficient for large spectral data sets, such as those produced by mapping 
spectrometers (NIMS,AVXBIS, etc. ) 
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Figure 2. Actual and calculated bi-directional reflectance spectra for the anorthite-enst&tite mixture series. The dotted 
curves, labelled AN for anorthite and EN for enstatite, are the endmember spectra used to model the mixture spectra. 
The smooth curves are the actual and calculated spectra plotted together. The number located in the 0.9 um. absorption 
band corresponds to the % mass of anorthite in the mixture. Each pair of calculated and actual spectra are offset for 
clarity and to emphasize the quality of the curve fits. At this scale there is almost no difference between the claculated 
and actual spectra. 
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Figure 2. Actual and calculated relative geometrici&l 
cross sections (f-parameters) of anorthite and enstatite 
plotted against one another for the mixture series shown 
in Fig. 1. The calculated values are obtained from the 
curve fitting technique described in the abstract. The 
calculated values are generally within 1.5% of the actual 
values which is adequate for estimating mineral 
abundances in the mixtures. 
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Figure 3 . This is the same type of plot as in Fig 2 for the 
olivine-magnetite mixture series. The letters correspond 
to the mass fraction ratios where A — 25/75, B = 50/50, 
Css 75/25, 90/10, E-95/5, and the first number is 

the % mass of olivine in the mixtures. The lines between 
data points were added to emphasise the systematic 
©verpredietkm of olivine content md underpr® diction ©f 
magnetite content. . 
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